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Abstract 

We develop a field-theoretical description of dynamical heterogeneities and fluctuations in supercooled liquids 
close to the (avoided) MCT singularity. Using quasi-equilibrium arguments we eliminate time from the description 
and we completely characterize fluctuations in the beta regime. We identify different sources of fluctuations and 
show that the most relevant ones are associated to variations of "self-induced disorder" in the initial condition 
of the dynamics. It follows that heterogeneites can be describes through a cubic field theory with an effective 
random field term. The phenomenon of perturbative dimensional reduction ensues, well known in random field 
problems, which implies an upper critical dimension of the theory equal to 8. We apply our theory to finite size 
scaling for mean-field systems and we test its prediction against numerical simulations. 

1 Introduction 

The heterogeneous character of glassy dynamics has been object of extensive study in the last decade [T] . Ex- 
periments, simulations and theory converge to a description of supercooled liquids where, on approaching the 
glass transition, relaxation requires cooperative motions on high mobility regions of increasing size and life time. 
An important theoretical step in the understanding of dynamical hetrogeneities has consisted in realize that the 
current theory of glassy dynamics, the Mode Coupling Theory (MCT) [2], predicts a growing dynamical length 
as the Mode Coupling critical point is approached. This was first noted in the context of disordered mean field 
systems where MCT is exact , and later confirmed with diagrammatic approaches to the dynamics of liquids |4] . 
In the resulting picture, the dynamical heterogeneities are captured by a time dependent four point correlation 
function, whose associated dynamical length diverges at the Mode Coupling critical point. As it is well known, 
this divergence, which is genuine in mean-field, is in real systems an artefactual consequence of MCT that neglects 
activated processes. The divergence is cut-off as the MCT dominated regime at high temperature crosses over 
to the barrier dominated regime at low temperature. With this caveat, the MCT prediction of a pseudo-critical 
growth of dynamical correlations has been largely confirmed in numerical simulations [5] and experiments fBl. 
However, corrections to MCT are at work as soon as the mean-field approximation is not exact. Two kinds of 
corrections to MCT can be expected: those due to critical fluctuations which arc not well described by mean-field 
theory, and those due to barrier jumping processes. Clarification of both kind of fluctuations is necessary to have 
an accomplished theory of glassy dynamics. Unfortunately both kind of phenomena are poorly understood. 

The goal of this paper is to present an in-depth analysis of perturbation theory around MCT to study critical 
fluctuations. In doing that we will neglect barrier jumping which is intrinsically of non-perturbative nature. 

The Mode Coupling (MC) approximation describes an ergodicity breaking transition where a system prepared 
in an equilibrium initial condition remains confined in its vicinity. Correspondingly, two point connected corre- 
lation functions develop an infinitely long plateau. This ergodicity breaking can be interpreted in the broader 
perspective of Random First Order Theory [7 . This theory predicts that within the approximations in which 
MCT is valid, at the dynamical transition the space of equilibrium configuration is partitioned in in an exponen- 
tially large number of metastable states. Several aspects of dynamical freezing can then be conveniently studied 



through equihbrium techniques, introducing appropriate constraints in the Bohzmann-Gibbs measure 8 • The 
free-energy as a function of the constraints provides a purely static field theoretical description of the MC er- 
godicity breaking transition. This description has indeed been crucial to the first theoretical recognition of the 
growth of a dynamical susceptibility at a MC transition [3]. In this paper, we exploit this constrained equilibrium 
technique to devise a theory of critical fiuctuations. The various dynamical characterization of fluctuations will 
be expressed in reparametrization invariant form eliminating the time dependence in favor of a dependence on 
the average value of the (two point) correlation function itself^ This perspective allows enormous simplification 
with respect to the dynamical perturbation theory |10| which at present is limited to the gaussian approximation. 
Previous studies have stressed the importance of emerging reparametrization invariance at large times as a soft 
mode of fluctuations in jlli glassy dynamics. Our approach will allow us to give a universal description of these 
modes in the beta regime where dynamical correlation functions are close to their plateau value. 

The main thesis of this paper is that reparametrization invariant fluctuations for temperatures close to the 
mode coupling critical temperature Td and values of the correlations close to the plateau value can be described 
in terms of a field theory of the kind 

H[<j>Mx)] = Jdx ^{\/^{x)f + ie + Se{xmx) + (1) 

where (j){x) is a local fiuctuation of the overlap away from the plateau value, e = T — is the deviation from 
the critical temperature, 5 is a coupling constant and Se(x) is an effective random temperature term, distributed 
with gaussian statistics and delta correlated in space. The effective Hamiltonian ([T]) coincides with the one that 
describes the spinodal point of the Random Field Ising model (RFIM) [12]. We find that both problems are 
perturbatively in the same universality class. The random temperature term is the ultimate consequence of 
dynamic heterogeneity and is a formal expression of "self-induced disorder" sometimes advocated to describe 
structural glasses. The role of this term is crucial. Random field models are well studied systems. It is well 
known that the random field changes the singular behavior of the theory. In particular in perturbation theory 
one finds the phenomenon of "dimensional reduction" which states that the singularities of the random model in 
dimension D are identical to the ones in absence of disorder in dimension D — 2. It follows that the upper critical 
dimension above which fiuctuations can be expected to have a Gaussian nature is found to be eight rather then 
six as it could expected from a pure 0"^ theory. It remains to find out if the Random Field Ising Model has a 
relevance for glassy dynamics beyond perturbation theory in the barrier dominated regime. 

The rest of the paper is organized as follows: in section[2]we analyze the sources of fiuctuations in the systems 
and we define correlation functions sensitive to them. In section |3] we discuss constrained measures. We explain 
their use in the computation of correlation functions and how to obtain them from replica field theory. In section 
|4]we analyze the replica field theory close to the MC critical temperature and study the quadratic fluctuations. 
In section [5] we analyze deeply perturbation theory and we derive the effective fleld theory ([l} . Section |6] is 
devoted to flnite size scaling in mean-field systems. The results of this last analysis are compared with numerical 
simulations in section [T] Finally we expose some concluding remarks in [51 

A partial account of the theory and simulations exposed in this paper has been given in [T31 . 

2 Measures of Fluctuations 

The theory exposed in this paper will be largely independent on the choice of systems. The main hypothesis we 
will make is that in some approximation a MC transition is present and we will study the generic behavior of 
fluctuations around it. Our theory apply equally well to describe critical fiuctuations around the avoided MC 
transition in liquids as well as finite size scaling around MCT in mean-field spin models where the transition is 
sharp in the thermodynamic limit. With the former application in mind, in the following we will use the language 
of field theory. In our formulae finite size scaling in Mean-Field models can be obtained simply replacing all space 
integrations by an overall volume factor N. 

For notational convenience we will represent the systems in terms of spin variables fixed in space Si = ±1, 
i = 1^...^N. With this notation we can equally well describe genuine spin systems like spin glasses, but also 
liquid systems in a lattice gas representation where we divide the volume in small cells and use the spin -taking 

^The terminology is mutated from asymptotic aging theory where time dependence is expressed through dependence on average 
correlations [9]. 
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the two values ±1- to represent the occupancy of the cells0 We will use as order parameter of freezing the 
correlation function, or overlap, among spin configurations. Given two configurations of the system S and S" we 
can define the local value of the overlap coarse-grained over some volumes v containing a large number of spins 
|u| S> 1, qx{S,S') = |w|~^X]iet,^ SiS'^. Different notions of correlations among configurations e.g. the one used 
in |15j lead to the same results, modulo a redefinition of the non- universal parameters appearing in ([Ij. If we 
denote by S{t) the configuration of the system at time t, the time dependent correlation function can be written 
as C{t,0) ~ y Jy dxqx{S{0), S{t)). The objects of our analysis will be the fiuctuations in the global quantity 
C{t,0) and the local quantities qx{S{0), S{t)) as they can be studied through 4-point or higher order correlation 
functions. 

We would like to separate the contributions of different source of fluctuations of C{t, 0). For structural glasses 
we would like to distinguish fluctuations among different trajectories that start from the same initial configuration 
from fluctuations due to changes in the initial condition itself. Recent numerical studies in supercooled liquids 
have emphasized the importance of this separation to study the influence of the structure in the development 
of dynamical heterogeneities |14) . For systems with quenched disorder, like e.g. spin glasses, one has a third 
source of fluctuations in the choice of the quenched interactions. In the following we assume without loss of 
generality the presence of some quenched disorder. If there is no disorder the respective averages are immaterial. 
We denote by (•) the average over trajectories that start from the same initial condition. This was called iso- 
configurational average in |14j . The iso-configurational average can be the average over the initial velocities in 
the case of Newtonian dynamics or the average over thermal noise along the trajectories in the case of stochastic 
heat bath dynamics. The initial condition is denoted by S{0) = S'^ and will always be chosen as an equilibrium 
configuration in this paper. The corresponding average will be denoted by |-]. Finally averages over quenched 
disorder will be denoted by E(-). A widely used measure of dynamical correlations is the 4-point correlation jl6] 
X4(i) — A^IE|(C(i)^)] — (E|(C(i))]) . In order to quantify the contribution of each source of noise to this function 
we use a decomposition of X4 in three different terms Xi — Xth + Xhet + Xdis defined as [17] 

^Xth{t) = E(I(C(t, 0)2)1) - m{C{t,Q)?l) 
^XHet{t) = E([(C(t, 0))^!) - E([(C(t, 0))f ) 

^Xa^s{t) = E(I(C(t, 0))f ) - E(I(C(t, 0))!)^ (2) 

These susceptibilities are the space integral of correlation functions that we will denote respectively Gth{x,t), 
Ghet{x,t) and Gdis{x,t). For example Ghet{x,t) can be expressed as: 

G;,et(x,t)=E(I(fe(5(0),5(t)))(go(^(0),5(i)))l-[(q,(5(0),5(t)))l[(go(5(0),5(0))D (3) 
In the case of liquids where quenched disorder is absent one has Xdis = and 

^xtH{t)^i{Git,or)j^i{Git,o)rj 
^xHetit)^i{cit,o)rj-i{Git,o))r. (4) 

In the following we will analyze the behavior of these three characterizations of fluctuations and predict their 
behavior for times such that the average correlation function Gav{t) = E(|(C(<, 0))]) is close to the plateau value 
Gp. As we will see in next section this can be achieved through quasi-equilibrium techniques at the price of 
eliminating time from the description. In the aforementioned time regime Gav (t) is a decreasing function of time. 
We can express time dependence through the dependence on Gavit) itself. For any time dependent quantity 0{t) 
we write 0{G) = 0{t)\c^^(t)=c- All the time dependence is condensed in the dependence of Gav{t) on time that 
we will leave unspecifled. 



^Having in mind a monodisperse systems occupying a D dimensional box of linear size L, we can divide the volume in N = (L/ a)^ 
cells of linear size a of the order of a fraction of the particle diameter. We then assign to each cell i the variable Si which takes the 
value 1 if the center of a particle lies in the box and the value -1 otherwise. 
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Figure 1: The behavior of Xth for the spherical p-spin model with p = 3. The curves are at temperatures 
T = 0.665, 0.68, 0.695, 0.71, 0.725, 0.74, 0.77, 0.785, 0.8 the critical temperature being Td = ^3/8 = 0.612. Upper 
panel xth as a function of time. Inset: C as a function of time. Lower panel xth as a function of C. 



3 Quasi-equilibrium in Dynamics and Constrained Boltzmann-Gibbs 
measures. 

In this section we discuss how to obtain information about equilibrium dynamics through the use of constrained 
equilibrium measures. This possibility relies in the phenomenon of time scale separation observed in glassy 
dynamics, where one can separate the degrees of freedom in fast and slow ones. 

The dynamics of liquids close to the glass transition can be described as a slow process where the system passes 
from one metastable state to another. Time scale separation tells us that approximate equilibrium establishes 
in a given metastable state before a new state can be found. The equilibration time within a metastable state 
is identified by the time that the correlation function takes to stabilize to the plateau value Cp. In the beta 
regime metastable states are sampled in a quasi-ergodic fashion. On this time scale, the different four point 
correlation functions introduced in the previous section can be then evaluated using constrained equilibrium 
measures that select the relevant metastable states. The set of constraints to be introduced should insure that 
the relevant regions of configuration space in the restricted measure coincide with the ones sampled by the 
dynamics. The simplest possibility is to impose that in each region of space the overlap with a well thermalized 
initial condition takes a fixed value. We will suppose that this specification of the local overlaps provides is a 
sufficient determination of the metastable states and assume that configurations that have a fixed overlap close 
to the plateau value with an equilibrium initial condition are sampled (almost) ergodically. This hypothesis 
-sometimes called separability [18]- can be checked directly in mean-field spin glass systems and we believe to be 
valid in supercooled liquids. In fact, we expect it to apply every time that glassiness is caused by the ruggedness 
of an energy landscapqfl. In the passage from dynamics to this quasi-equilibrium description we loose of course 
the possibility of studying the time dependence of the various quantities, that, as mentioned in the previous 
section, will be expressed instead as functions of the overlap in a time reparametrization invariant representation. 
In figure [1] we illustrate how the four point dynamical susceptibility of the spherical p-spin model [3] looks like if 
we operate this change of perspective. 

Let us remark at this point that the use of a time independent description of dynamical quantities has been 

''On the contrary, we do not expect to apply in systems like kinetically constrained model, where the Hamiltonian is trivial. In 
this case the overlap does not give a sufficient determination of the metastable states |19| . 
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widely used in the theory of aging [5] , where reparameterization invariance emerges as an asymptotic continuous 
symmetry at large times. It was then proposed that this asymptotic zero mode could be used to characterize 
fluctuations 11 in glassy states. Being based on symmetry considerations this theory is very general, and concerns 
features of both the beta and the alpha regimes. Our theory, being based on quasi-equilibrium considerations is 
less general and more specific: it will enable to give a detailed description of the beta regime, but it is limited to 
that. We will see however, when discussing simulations, that looking at the data in reparametrization invariant 
form is useful and inspiring also in the alpha regime. 

In the rest of the paper we will concentrate on values of the temperature close to T^, and ignore the possibility 
of an ideal glass transition at a lower temperature Tk- 

We will concentrate on systems where either quenched disorder is absent, like in real liquids, or, if disorder 
is present its effect is weak and physical quantities can be evaluated to the leading order by the "annealed 
approximation" . This is a stronger property than the usual self-averaging property of the free-energy and states 
that the partition function has small sample to sample fluctuations. Systems of this kind are often used to model 
structural glasses, and include, among others, fully connected p-spin models, spin models on diluted random 
graphs and finite range mean-field models in the Kac limit. 



3.1 Effective Potential: a Landau field theoretical functional for the glass transition 

According to the discussion of the previous section we can select metastable states just choosing random equi- 
librium configurations S'^ and restricting the Boltzmann measure to configurations that do not differ too much 
from 5°. We achieve this fixing the local overlaps qx{S, S*") to preassigned values px and defining a constrained 
measure [5] 

For systems that are separable in the sense specified in the previous section, the equilibrium metastable states are 
selected choosing in all points of space Px — Cp, profiles that deviate from this shape allow to probe fluctuations. 
The partition function Z[S^,Px] is in fact directly related to the probability of the overlap profile 

where Z is the unconstrained partition function. The large deviation functional W{px, 5*") depends on the choice 
of the overlap profile but also on the choice of the reference configuration S*" and on quenched noise in the case 
of disordered systems. In our formalism any dependence on 5'° quantifies the notion of "self-generated disorder" 
often advocated in the physics of structural glasses |2D] . Previous studies have concentrated on the average value 
oiW 8 . Depending on the nature of system under study, one can expect that the fluctuations of W with respect 
to S'° and J are more or less strong. For example, in a fully connected model is a function of a single global 
overlap parameter and self-averaging in the thermodynamic limit. Fluctuations decrease as powers of the system 
size. We will see however that fluctuations of the correlation functions and fluctuations in the potential can be 
related to each other. The entire probability distribution of W is therefore relevant to a complete description of 
glassy systems. 

In fact, the present formalism allows in principle to compute the dynamic correlation functions that we have 
defined in the previous section in reparametrization invariant form. To this scope we introduce the generating 
function of the overlap T(hx\SQ) defined by 

l3r(h^,S") ^ I -n-n p-PW{p^,S'')+f dxh^p. 



and define the static analogue of the correlations ([3]) in presence of the field hx as 

Gthix -y,h)= EiliPxPy) - {Px){Py)l) = E{l{pxPy)cl) 

GU^~-y,h)^Eil{px){Pyn - l{Px)E{Pym^m{Px){Py)ic) 

G*s(a;,/i) =E(I(p,)M(Py>l)-IE(I(p.>l)IE(I(p,)l) =E(I(p.)lI(p,>l)c (8) 
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where we have denoted here by (•) the equihbrium average in presence of and with a subscript "c" subtraction 
of the disconnected part. It is easy to check that the various correlations are related to the derivatives of the 
moments of the F functional according to 



Gthix - y,h) ^ 
Ghetix, h) = E 
Gdis{x, h) = E 



E 



S^Tih\S") 



ShxShy 



(Sr(/i|50) ST{h\s'>) 



Shx 

sr{h\s°) 



Shx 



Sh„ 



5T{h\S°) 



Sh„ 



(9) 



If we fix the field hx in such a way that E(|(p2:)]) — q in all points of space we get the correlation functions as 
a function of q. We need then a method to compute the the cumulants of the functional W or equivalently the 
ones of r. 



3.2 Effective potential and replicas 

It is interesting to compute both the average of the potential W(qx, S^) and its fluctuations. The replica method 
gives us a simple framework to undertake this task. As discussed many times jS], the average W^^\px) = 
MlW{px,S°)} can be computed considering 



(10) 



valid for integer mEI Notice that here the total number of replicas, which includes the reference configuration 
Sq and the m copies of the constrained system, is n = m + 1 . The free-energy functional is obtained from an 
analytic continuation to m = 0, i.e. the total number of replicas n tends to 1. 



W^'\px)^-T^^^ 
dm 



F 



(11) 



m=0 



where F is the average unconstrained free-energy of the system. Similarly one can get the second cumulants 



wtliiPx,p'x) = mwiPx,s'>)wip',,s'>) 



wj,2{Px,p'x) = mw{px,s'mwip^,s') 



a2iogE([z[50,pj"iz[50,p;]"i) 



dnidn2 

,712—0 

dHogE{iz[s",pxrqiz[s",p',r-}) 



dnidn2 



(12) 
, (13) 



ni ,n2— 



where in the second equation we have exchanged the logarithm and the average over the disorder, thanks to the 
annealed approximation. Higher order cumulants can be analogously obtained through more involved analytic 
continuations. 

In order to unify the notation and treat all cases in parallel it is convenient at this point to introduce the 
(formal) replica action S[Qx] for n replicas for fixed values of their mutual overlap Qab{x) (a, h — 1, n) from: 



= -E y 



n X{ma.l,{x)-qx{S'',S'')) , 



(14) 



{S"}S=i 



a,b=e X 



from which, integrating over some of the elements of the replica matrix and fixing some others one can get the 
moments of W . For example one has that W'^^^px) can be computed by a replica action with n — m+l replicas 
for m — > 0. Renumbering the replicas in a way that a = 0, 1, m one has 



/m 
I?gab(x)e-^['3==1 n S{QoAx) - P{x)) , 
o=l 



(15) 



■^Thanks to the hypothesis of self-averageness of the partition function 1/Z Ri 1/E(Z) the average over disorder in mUl l does not 
require additional care. 
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Analogously the correlation functions can be computed from a replica action with respectively n = ni + 712 + 1 
and n = rii + ^2 + 2 replicas for ni, 7^2 — > 0. 



d 



dnidn2 
d 



2?Qa6(x)e-^[Q^l Y[ SiQoAx) - p{x)) n ^{Q^Ax) - P'i^m) 

-iv-^-- a=l a=ni + l 

/ni ni+n2 
2?Qa6(x)e-^[«^l Y[ SiQoJx) - p{x)) \{ 5{Qo.Ax) - 
-iv-^-- a=l a=rii + l 



dnidn2 



where in the first case we have renumbered the replicas in a way that a — 0,1, rii + 712 and in the second 
a = 0, 0', 1, ni +712. 

We would like at this point to remind that in disordered mean-field models there is a close relation between 
the Mode Coupling dynamical transition and the shape of W^^'>{q), which in that case is a function of a single 
variable. In fact, the transition temperature Td looks as a spinodal temperature for the potential VF^^^ This has 
a single minimum at low values of q at high temperatures, and develops a second minimum right at Td for the 
value of the overlap q — Cp. 

We argue that in a separable system, where the measure ^ correctly samples metastable states, this is 
the generic situation. If metastability is found in some dynamical approximation, an approximation for statics 
with the same physical content should lead to the appearance of a secondary minimum in the average effective 
potential corresponding to the constant profile Px — Cp. Recent analysis of MCT as a Landau expansion [21] on 
one hand and reproduction of MCT results from replica Orstein-Zernike equations |22j on the other corroborate 
this point of view. 



4 The replica action close to 

We enter now in the core of our analysis, and we study fiuctuations for theories S{Q) that at the level of 
homogeneous {i.e. space independent) saddle point exhibit a dynamical phase transition at a temperature Td. 
This is associated to the appearance of a horizontal inflection point at Cp in the effective potential W'^^^ (q), which 
becomes a minimum below Td. As explained in detail in [5] (see also |23| ) this inflection point is described by a 
n = 1 replica symmetric saddle point where Qab{x) = Qfi, = Cp for all x and a^h. We wish to describe overlap 
fluctuations for T in the vicinity of Td and Px in the vicinity of Cp. A natural point of expansion of the action 
S[Q] is the homogeneous saddle point just described. We can then expand in e = T — and the difference of 
Qab(x) with Cp, 4)ab{x) — Qab(x) — Cp loT a ^ b and ipaa = 0. To the leading order one has a cubic theory 

ab ab 

+ j dx dy (l)ab{x)Mat-cd{x,y)(j)cd{y) 

ab\cd 

+ j dx dy dz ^ ^ab;cd;ef{x,y,z)(j)ab{x)4>cd{y)4'ef{z). (18) 
ab\cd\ef 

The second term vanishes for n = 1 where Q'^ is a solution to the saddle point equations. For generic n however 
this is not the case, this will be a term of order n — 1 that has to be kept in the expansion. To the lowest order 
in a gradient expansion and rescaling the variables to reabsorb superfluous constants, the action reads 

S[Q] = SolQ"] + j dx Y.{AaM + e)<Pab{x) 

ab 

+ ^(^ V(/ia6(a;))^ ^ ^ ^ 't^'^b[x)Mab-cd4'cd[x) 
ab ab\cd 

■cd;ef(l^ab{x)(l)cd{x)(j)ef{x). (19) 

a&;c(i;e/ 
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Notice that the temperature couples hnearly with (j)ab{x). This is due to the choice of the point of expansion as the 
inflection point at Td- The components of Aab^ the mass operator Mab-cd and of the bare vertex Tab;cd;ef reflect 
the symmetry of the saddle point and should then depend only on the number indexes that are equal or different. 
This immediately imply that the coefficients Aab for a ^ b should then be all equal Aab = A{n) '-^ (n — 1)A, and 
that the quadratic form can be written as 



\ ab ab a b ab 

The inclusion of all possible replica symmetric cubic vertexes gives rise to a cubic part p4] : 

SsiM = Jdx £(3) 

OJl ^ 4'ab(f>bc(f>ca + ^2 ^ (f>lb + 



(20) 



(21) 



£(3) 



ibc 



+ 4'lb(f'ac + ^^4 ^ 4'lb'^>cd + ^5 X! '^ab(t>ac't'bd + 



ibc 



abed 



abed 



ibcd 



abcde 



■^8 ^ 4>ab(t^cd4>ef 
abedef 



however, we will show that only the first two terms are relevant for n — > 1. 

Notice that the average potential within mean- field is evaluated by a saddle point (pabi^) 
inserted in (fTO]) gives, to the lowest order in n — 1, 



(22) 



(t>{x), which, 



g = ^(^2 - wi) 

6 



dx i(V0)2 



1 

—mi c 
2 



(23) 



Since by hypothesis we have developed around the horizontal inflection point for e = 0, we must have mi — 0. 

Of course, different choices for the point of expansion are possible. In the region T < it is convenient to 
expand around the replica symmetric saddle point Qab = Qea{T) at temperature T, i.e. around the point that 
describe the minimum of the average effective potential (|23p . This choice leads to an action like (|19l) where all 
terms linear in are absent. The factor mi is in this case non zero and proportional to \f~k. With this choice 
of the expansion point the average potential reads VF^^-* ((/ix) = Vl^(O) + J dx i(V(/))^ + ^micfP + gcj)^. In the 
following we will use both expansions, without introducing separate notations for the two. 

Notice that we have written the expansion (|19p for generic D dimensional extended systems. However, as 
we will see in section ([6]), the same expansion can be used to describe finite size corrections in fully connected 
disordered models and models on diluted random graphs. In that case the various overlaps are global quantities, 
the gradient term is absent and space integration is just substituted by an overall multiplication by the system 
volume N. 



4.1 Quadratic free-energy fluctuations 

In this section we discuss the correlation functions at the quadratic (one loop) level, neglecting the cubic part 
of the action. We study the fluctuations of the effective potential with respect to the choice of the initial 

(2) (2) 

configuration and the choice of quenched disorder. Our task is to compute Wf^J^ and W^^^. To this scope we 
observe that the one loop order can be estimated as the saddle point of the ni, n2 derivative of S'[0] where some 
of the matrix elements are fixed. Let us start from Wf^J^{(j), (j)'). In this case we need to consider a saddle point 
of the action with n = 1 + ni + n2 replicas where the elements 0o,a are fixed to 4>o.a{x) = (l){x) for a = 1, ni, 
4'o,a{x) = (f)' {x) for a = ni + l,...,rii + 712. In presence of such constraints it is natural to look to a saddle 
point which is symmetric with respect to all the permutations that leave invariant the values of constraints, i.e. 
the group S'„j x of independent permutations of the first group and the second group of replicas among 
themselves. This is parametrized in terms of three fields i^{x), tp'{x) and ^^0(2;) such that all the couples of 
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replicas a,b — 1, ni have the same overlap (t>ab{x) = ^{x), all the couples of replicas a,b = ni + \, ni + n2 
have the same overlap <f>ab{x) — ip'{x), and all the couples of replicas a — l,...,ni b — ni + l,...,ri2 have an 
overlap 4>abix) = 4'o{x). Inserting this ansatz, one realizes that in the leading order in ni,n2 — the equations 
for ■0 (resp. ip') are independent from (/)' and "00 (resp. (f) and V'o) and coincide with the ones that appear in the 
computation of VK*-^-*. In the limit of small e the solution is simply — (/), ip' — (/)' and -00 = 0, which gives 

W^/.ei[<^-<^^] = IW'{0\S")}+A J dx (0(a;) +0'(x)) - {m^+m^) j dx 4>{x)<f>' (x). (24) 

This formula has a clear interpretation: the effect of the heterogeneity in the reference configuration S*^ can 
be parametrized in terms of a space dependent random free-energy shift a{x) and a random temperature Se{x) 
which couples linearly to (j). This suggests that the potential W{(j)\S'^) can be written as 



W{<p\S°) = / dx 



-(V<?!))2 + (e + Se{x))<j>(x) + g(l){xf + a{x) 



(25) 



The free-energy shift and the random temperature are Gaussian mutually correlated fields with 

la{x)a{v)l = lW\0\S''mx-y)/V (26) 
I5e{x)6e{y)\ = -{m2 + m:i)5{x - y) (27) 
la{x)5e{y)l = A6{x-y) (28) 

where the consistency of the theory requires m2 -I- ni^ < 0. Formula (j25p is the central result of our paper, 
derived here at the level of quadratic free-energy fluctuation. The effective field theory for the dynamic glass 
transition coincides with the one describing the spinodal point of the Random Field Ising model and therefore 
both problems are in the same universality class. It is well known that the leading singularities of random field 
theories in perturbation theory are given by the tree diagrams [25j , or by the formal solution of the stochastic 
differential equation 

- Act) + {e + Se{x)) + 3g(j){x)'^ = 0. (29) 

In fact, even when this equation does not admit real solution, the complex solutions gives rise to a perturbation 
series for physical quantities with real coefficients. Though the analysis of the quadratic fluctuations give a strong 
hint about the validity of ([25|) . one could doubt that the inclusion of the vertexes in the theory modifies this 
result. As the matter of fact, in the next section we will analyze in depth the perturbation theory for T < Td 
and confirm (1251) to all orders in perturbation theory. 

(2) 

Let us now briefly turn our attention to W^^^. In order to compute this quantity one may follow a route 

(2) 

similar to the computation of Wf^J^. In this case however one can note that in annealed models the replicas 
and 0' have zero overlap, and an expansion around q = Cp is not justified. Rather, one has that the mutual 
overlap between the replicas labeled 1, rii and the ones label- led rii + l, ...,ni + n2 should be put to zero. This 

(2) 

leads to decoupling between the two groups of replicas and W^-^ = to the leading order. This is a remarkable 
result, showing that fluctuations due to disorder are much less important than fluctuations due to "self-generated 

(2) 

disorder", seen here as heterogeneities in the reference conflguration. We stress that the vanishing of W^^^^ 
at the quadratic level is consequence of the annealing hypothesis and certainly would not be true in systems 
where disorder fluctuations in the partition function have to be taken into account. In our view this absence of 
dependence on quenched disorder strongly supports the validity of long-range p-spin and similar models as good 
mean-fleld models for the structural glass transition. 

Let us now exploit ([M)) to compute the overlap correlation functions ((S]). First of all we notice that at the 
tree level calculation the potential W and the generator F are related by 

ri[h] = maxVF(i)(g) - J dx Kq,^. (30) 



This imply that the correlation function Gth at the one loop level is given by 

5Ti[h] _ 

ShxShy \ 5qx5qy 



Gth{x ^ y) = —— ^ { ——— \ . (31) 
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We also notice that to the same accuracy 

^^Hlt[h.h'] = wi%,q') (32) 
where q and q' are the maximizers of eq. pop for field h and h' respectively. A simple computation of the 

(2) 

derivative of F^^j shows that 

^''-(^ -y^= ^-^B^ = I '^y' - x')G,,{y ~ y') (33) 

oh^dh'y J oqx'dqy, 

(2) (2) 

The same computation for Gdis{x — y) would yield the same formula with W^^j at the place of Wf^J^, but as we 
have remarked Vy^jj = to the leading order. This implies that no singularity of Gdis can be detected in the 
quadratic theory. Beyond the gaussian approximation, the singularity of Gdis , if any, should be weaker then the 
one of Gth and Ghet and the possible presence of quenched disorder does not affect the universality class of the 

system. Eq. p3p shows that as soon as i-jf^jj^-^) is non zero, the order of the singularity in Ghet{x — y) is the 



double of the one of G 



th(x - y). Notice that ^^7^^^^ is precisely equal to Pe^^Jej^] = -(7712 + m3)6{x - y). 

We find the announced result that the largest source of fluctuations in the system comes from the heterogeneities 
in the initial condition. Its effect at the one-loop level is to double the singularity due to thermal fluctuations. 
If we specify to the form (1^^ we find 

°-^—M = 5[x-y){-^ + &g<l>{x)). (34) 
oqxoqy 

Fixing now a constant overlap profile in space we find that in momentum space 

Gth{kA) = \ ^ (35) 



n (h A\ (ma -l-ma) 

= -(P + 6#)(P + 6#)- 

Both propagators are singular at = and fc = 0, this corresponds to the divergence of the fluctuations at e = 
and q = Cp. The fluctuations for e 7^ can be obtains inserting for (f> a cut-off value of the order of the plateau 
(I^EA ^ y/\e\/g and 

Gthik) = ^ (37) 

G,^,ik) = ^ . (38) 

Notice that, within the present Gaussian approximation, the intensity of critical temperature fluctuations 
l^ejj^Cy] = — (m2 -I- m3)6(x — y) can be measured through the p ratio 



^ Ghetjk, 4>) 

~ Gth{kA) 



2 



(39) 



at the critical point e = 0, = 0. Beyond the Gaussian approximation, while and (1551) provide a clear 
indication that Ghet is more singular than Gth, it is not clear to us if the relation Ghet ~ continues to hold. 
This different scaling is an important result of our theory. Though our derivation is restricted to the beta regime, 
we will see in section ([T]) that numerical simulations show that a different scaling is also observed in the alpha 
regime. 



5 Perturbation Theory 

We would like to confirm the description of fluctuations through the potential ([^5)) by perturbation theory. In 
order to have a well defined point of expansion in perturbation theory, we assume T < Td and we expand around 
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the replica symmetric minimum of the action Qab{x) — qEA{T) for ah a, h. Our starting point is then 

\ ah ah ah ah J 

with £(3) [^] 

given by (j22p and toi ^ V— e. 
5.1 The bare propagators 

Let us now reobtain the results on the 4-point functions of section (|4.ip by studying the bare propagators of 
the replica field theory, as first derived in |26| . The analysis of a generic mass matrix with replica symmetric 
structure has been performed long ago by De Almeida and Thouless p7]. To analyze our case, we need to 
transpose their results to the case n — >■ 1. There are three in general distinct eigenvalues of the quadratic form 
named longitudinal (L), replicon (R) and anomalous (A) in the current terminology. For future reference we 
remind that the longitudinal sector correspond to fiuctuations such that 0^;^ = (/> independent of a, & (a 7^ 6), the 
replicon sector correspond to fiuctuation such that is vanishing for all a and the anomalous sector is the 

linear space of fiuctuations orthogonal to the previous two. In terms of the parameters mi, TO2 and ms of the 
quadratic form in (j40p . the eigenvalues for generic n read 

Afl(fc) = fc^+TOi (41) 
AL(fc) = fc^ +mi + (n- l)(m2 +nm3) (42) 

A^(fc) = fc2+mi + i^i^TO2 (43) 

The replicon Ai(;(0) — nii ^ -\/— e is critical at the transition. Notice that for n = 1 the longitudinal eigenvalue 
becomes degenerate with the replicon and therefore it also become critical at the transition. 

It is already clear at this point that the replicon and longitudinal sections will give the most singular con- 
tribution to the perturbation theory. The propagator matrix of the theory in momentum space, G^ah-cJ^) ~ 
(A:^ + M) ^.^^ has the same replica symmetric structure as the mass matrix, so that we can write 

^(0) n X ifiachd + ^adhcl , {Sac + Sad + he + 5bd) , , , 

Gab,cdik) = 91 ^ + 92 ^ + 93 (44) 

The coefficients gi , 52 and (73 can be easily expressed in terms of the eigenvalues of the mass matrix, or in terms 
of the m parameters, the result for n = 1 is: 

51 = rjT (45) 

fc^ + mi 

g, = I ^— (46) 

m2/2 + k'^ + rni k? + mi k"^ + toi 



1 



93. 



fc^ + mi 



-(m2 + m3) m2 



fc^ + TOl mill + /c^ + TOl 



(47) 



The replica formalism naturally embeds the distinction among different sources of fiuctuations discussed in section 
^ and allows to easily compute the propagators Gth and Ghet of section [5] to the one loop order. This can be 
done noticing that E(|((^^)]) = Gab,ab, E(|((^)^]) = Gab,ac where all the replica indexes are different one from 
another 

Gth{k) = Gab-ab — Gab-ac = 91 — g2 ^ 7^5— (48) 

fc^ + mi 

n (U\ n ^ -(mz + m3) , - 

Ghet(k) = Gab-.ac = g2 + g3 ^ , , ■ (49) 

(fc-^ + TOi)^ 

Coherently with the results of section HTTl the singularity of Gth is a single pole while the one of Ghet is a double 
pole. Within the replica formalism the origin of the double singularity stems from the degeneracy of the replicon 
and the longitudinal eigenvalues for n — ^ 1, which both become critical at the transition. 
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5.2 A perturbative derivation of the Stochastic Equation 

We would like now to analyze the complete theory and show that the leading singularities in perturbation theory 
coincide with the one given by the stochastic differential equation (j29[) . The analysis we perform is similar 
to the analysis of random field models as originally put forward by Parisi and Sourlas [25] . This is based on 
dimensional evaluation of the various vertexes of the theory for n — > 1 exploiting a change of basis that generalizes 
the one suggested by Cardy in [28^ for the RFIM. We note first of all that the leading singularities come from 
the replicon and longitudinal modes, that become critical at the transition. We therefore concentrate on these 
modes, for which the matrix (jjat is such that is independent of a. Symmetric matrices with this property 

form a linear space of dimension (n — l)(n — 2)/2. We now describe fluctuation in a different basis separating 
the fluctuating replica matrix in replica symmetric part, independent of the indexes, plus a replica symmetry 
breaking fiuctuations and write, for a ^ b and all points in space, 

0ab = - + Uab ^ + Xah (50) 

^Xab = Va (51) 

h 

^XahUab=^ (52) 
a, 6 

Where we have defined IJab as a constant block matrix which has all elements equal to zero except the ones for 
which & = a + 1 and a is odd or 6 = a — 1 and a is even, i.e. 

\ otherwise ^ ' 

In this new basis we will be able, on the one hand, to perform the n — > 1 limit directly in the action and, on the 
other hand, to evaluate the scaling dimension of the different terms in the action in order to keep only the most 
singular ones. Notice that the "vectors" <j>ab defined by ([52]) span the longitudinal and replicon sector and that 
the matrices Xab verifying the above relations form a linear space of dimension (n — l)(n — 2)/2 — 2. We observe 
now that 

^cl)at = ^ + in-l){(j)~^uj) (54) 

b 

€b = n2c.0 + Y.^lb + n(n - - |)2. (55) 

a,b a.b 

Neglecting all the terms that vanish in the n — )• 1 limit, the quadratic part of the action reads 



52 = j dx 2(V0)(Vw)+^(Vxab)'+mi(</)w + i^Xa6) + ^(™2+m3)w' 

\ ah a.b I 



(56) 



Let us now study the mi-mass dimensions D^, D^^ and of the fields a; and Xab- As usual we impose that all 
terms in 5*2 have the same dimension. We consider the case in which m2 and 7713 remain finite at the transition, 
while mi — 0. In this case we can write 



D^ = D^ + l (57) 
1 
2 



D^=D^ + - (58) 
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Let us now analyze the vertexes. Expressed in the new variables the first four vertexes of (j22[) read: 
5^ (l^abMca = (n - l)(n - 2)(0 - |)3 - 30^ Xab + 3n(n - 2)u:<i>^ + Tr 

a,b,c ah 

+wTrx'C/ -\^Y. + 3(/.c^2n(n _ 2) + ~ 2) (59) 

ah 

J2 <Plb - n(n - 1)(0 - |)3 + 30 5] x^b + 3nc.(/.2 ^ ^ ^3^ 

+3^E^'fcC^ab-|^E^'fc + 6na.20 + n^ (60) 



(61) 



^<^2^0ac= (w + (n-l)(</>- it^)) n2w0 + ^ Xab + n(n - 1)((/)- 

a,b,c \ a,b 

X! <^Qb<^c<i = X! '^ab'/'ac (62) 

all the other combinations that appear in ^3, i.e. vertexes 5 to 8, give just rise to terms proportional to cj"^ plus 
terms that vanish for n ^ 1. Notice that all the terms proportional to cfr' vanishes for n — > 1. In the new basis, 
the mass dimensions of the vertexes that survive for n — )• 1 are 

(t)^xib 



b 

3 



XI ^a6> Tr X 

ah 



^3 



-> 3D^- 


■f 1 


(63) 




3 

^2 


(64) 




-f 2 


(65) 




-1-3 


(66) 



The leading singular behavior in perturbation theory of the theory for mi — > is dictated by the the first two 
vertexes (j63p . which are the ones of lower dimension. Neglecting therefore the subleading vertexes we find that 
we can write the action as 

dx -(m2 + ms)^^ +0; (-A(/) + micj) + "igcj)^) + - ^ Xab (-A + mi + 6g(l))xab (67) 

a,b 

where g = 002— u!i- We observe that the matrix field Xab{x) has (n— l)(n— 2)/2 — 2 — > —2 independent components 
and becomes equivalent to a couple of anticommuting fermion fields xi^) ^-nd xi^) for n — >■ 1. Equivalently, we 
can observe that the explicitly integration over the Xab fields gives rise to 

det (-A + mi + 6.g0(x))^"'""^^^""^^/'* — > det (-A + mi + 6g(f){x)) . (68) 

n— 

We finally recognize in the action ()67p a Parisi-Sourlas supersymmetric theory associated with the stochastic 
equation 

- A0(x) + mi(j){x) + 3g(j){x)^ + 5e(x) = (69) 

where de{x) is a gaussian field with variance 

lSe{x)Se{y)l = -(ma + m3)S{x - y). (70) 

If we impose that all the terms in the action have the same scaling dimension we find = 1 and in D spatial 
dimension the action has dimension — -^+4. It is well known that Parisi-Sourlas actions present a supersymmetry 
that leads to the phenomenon of dimensional reduction. The perturbation theory of the system in a random 
field in dimension D coincides with the one of a pure system in dimension D — 2. Coherently the upper critical 
dimension of the theory is promoted to 8 from the value 6 of the pure 0"^ theory. 
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From the (formal) solution of eq. (I69p we can obtain the correlation functions Gth and Ghet through 

linear response theory, using (j70p 



Gtkix -y)^ = ^ Uix) X 6eiy)l 

de(y) TO2 + ma 

Ghetix~y)^mx)mic- (71) 

The derivation leading to ([ST]). ([55]) . ([701) is vahd within perturbation theory. Its expression should be 

considered as a formal writing valid within the perturbative context. In fact, due to the cubic nature of the 
potential, eq. (|69p admits a real solution only if Se{x) is sufficiently small in absolute value in all points of space. 
For a gaussian field Se{x) this condition is violated with probability one in an infinite space. The consequence of 
that is that the perturbation series should be divergent. In fact, the perturbation theory is formally identical to 
the one of branched polymers considered in 29 , with the crucial difference that in polymer's case the coupling 
constant g is purely imaginary |30j . while in the present case it is real. In the branched polymer case the 
perturbative series is resummable thanks to the fact that its terms have alternating signs. Here, all terms have 
the same sign, the resulting series is badly divergent and it seems hardly resummable. Recent work has used 
"Exact Renormalization Group" [31] methods to compute the exponents of the thermodynamic transition of the 
RFIM [32]. It is not clear to us if these methods can be useful to study the spinodal point. 

A fast way of realizing that perturbation theory should be divergent comes from considering equation (|69p 
in the homogeneous limit of space independent quantities. As we will describe in the next section, this allows to 
describe finite size scaling in mean-field models. In this case the equation reads 

mi(/) + 3g(/)2 + 5e = (72) 

where now the variance of Se is given by 

ISeSel^J-^^^I^. (73) 



The correlation function Ghet can be formally evaluated from the solution cj)* = i2g(5e 

Ghet{e,N) = ^{Mn~m') (74) 

This correlation can be evaluated in an expansion in Se, giving rise to well defined series with real coefficients. 
However it is clear that the series cannot be convergent as the averages in (|74|) receive an imaginary contribution 
from the square root in the region Se > m\/12g. 

Putting these problems of convergence aside, one can ask if, for typical disorder realization Se{x), the stochastic 
equation ([M)) could be used beyond perturbation theory in a description of the barrier jumping processes in a 
reparametrization invariant way. While we do not have a definite answer in general we will discuss the consequence 
of this hypothesis for mean-field finite size scaling in the next section. 



6 Finite size scaling around mean field 
6.1 The beta region 

The theory that we have developed can be easily generalized to describe finite size scaling in mean-field systems 
which have a genuine MC dynamical phase transition in the thermodynamic limit. These include disordered spin 
models like p-spin and Potts defined on fully connected or finitely connected random graphs, and problems that 
appear in computer science and information theory like the K-SAT or XOR-SAT problem or error correcting 
codes. For large but finite sizes N the dynamical transition is cut-off. Finite size scaling should describe the 
cross-over to criticality ioi N co and e — ?► 0. 

A phenomenological theory of dynamic finite size scaling for disordered mean field models has been first 
proposed in |33j . In order to interpret numerical results, in Ref. j33j it was assumed a sample dependent critical 
temperature and the Harris criterion was used to derive scaling variables and exponents. Our results, using a 
fundamental theoretical description, confirm and rationalize that analysis as far as reparametrization invariant 
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quantities are concerned. As discussed above, the origin of the random temperature term is in our theory a 
consequence of dynamical heterogeneity rather than of quenched disorder. 

To analyze mean-field finite size scaling we follow the lines drawn in the previous sections, except that in this 
case overlaps do not display space, that is we use observables integrated over the whole system. The relevant 
replica action is identical to the one discussed in section 2] without the gradient term and with space integration 
substituted by an overall volume factor N. 

Repeating the analysis of the replica action that led to ([^5)) and in this case, we get a description in 
terms of a single variable effective potential W((f>\So) describing the total overlap fiuctuations around the plateau 
value, (j) ^ q ~ Cp 

W{(f>\S°) = ^(015°) + N{[e + Se]<j> + g<j>^ + a) (75) 

where in this case both Se and a are gaussian covaring variables of order 1/^/N. In terms of the parameter of the 
replica field theory, g = i(cj2-wi), |aa] = ;^|iy2(0|S'o)l, faSej = j^A, {Se^ = - jj {m2 + ms) . Correspondingly 
we get the Parisi-Sourlas action 

S = N [^I^^±^u^ + u{e + 3#2) + x6#x) ■ (76) 

From formula (|76p the properties of finite size scaling readily follow from dimensional analysis. We are interested 
to the behavior of the various observables in the critical cross-over region for — > oo and e — > 0. This is the 
region of variables such that all terms in the action are of order one, namely (j) ^ \/[ef, uj ^ e ^ iV~^/^ and 
X, X ~ N~'^^^ . This allow to identify as scaling variables x = ipN^/'^ and y = eN^^'^ in quantities that depend on 
size, temperature and overlap. Mutating the results of section |4T] on the behavior of quadratic fluctuations, we 
find that, for — and N — )■ cxs, xth ~ —7= and Xhet ~ ri- Trading N for e, this implies that in the cross-over 

region 

Xth{<l^, e, N) = N^'^f,^{<i,N''\ e7Vi/2) (77) 
Xheti^, e, N) = N^'^fHet{4>N^'\ eN^/^). (78) 

In order to match the singularities for N — > 00, the scaling functions fth{x, y) and fhet{x, y) verify fth{0, 2/) ~ 

fth{x, 0) ^ ^ and fhet{0, 2/) ~ ^1 fhet{x, 0) ^ ^ for a; — >■ -|-oo or y — -l-oo, while a finite limit should be expected 
at small values of the two arguments. 

It is interesting to study the behavior of the ratio p = Xhet/ Xth scaling window. In the N ^ 00 limit 

this ratio can be related to the variance of the random field by the relation: lim0_>.o liiii7v->oo p(e — 0, (j), N) — 
— (m2 + 7713) where the limits should be taken in the order. In the scaling window on the other hand p(e, 0, N) — 
fhet{x,y)/ fthix.y)'^. Consistency requires that limj^^oo //iet(0, 2/)//t/i(0, y)^ = -(mi -1-7712). As we have observed 
in the previous section the perturbative series of the scaling functions are badly divergent and can hardly be 
computed analytically. Thus, this is as far as we can go from the perturbative analysis of (j76p . Numerical 
verification of the scaling forms ((78|) will be the object of next section. 



6.2 The alpha regime: some conjectures 

In the deep alpha regime, where < is of order one in absolute value, quasi-equilibrium in the sense we were 
using so far does not hold, the configuration space with average correlation C < Cp is not sampled ergodically and 
our theory does not apply. However, this does not make it less interesting to look at the time reparameterization 
invariant part of the fiuctuations and on the contrary calls for new theoretical ideas to be developed. One can 
hope that finite size scaling of mean field models is a simple enough setting to put forward hypotheses to be 
tested in the general case. This section will be then by nature much more conjectural and qualitative than the 
previous ones. 

In order to understand finite size scaling in the alpha region we can conjecture continuity with the behavior 
in the scaling regime (j) ~ N~^^^ and quote the results of the analysis of dynamic Gaussian fiuctuation theory 
developed in [10]. This predicts that the peak xt of the dynamic susceptibility as a function of time, lies deep 
in the alpha region and scales as If we assume that this scaling holds in the whole regime C < Cp, we can 

^In IIUI this behavior was found to hold only for conservative dynamics. Later analysis showed that the same behavior also holds 
generically even in absence of dynamically conserved quantities | 34| . 
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match the form (|78p if we suppose that the function fhet{x, y) behaves quadratically for large negative values of 
X. In this way Xhet{(f', N) = N^'^ fu^MN^'^ ,eN'^'^) ^ N(fighet{eN'^'^) for -(pN^/^ » 1, where the function 
ghet{y) should behaves as y^^ for large argument and take finite value for y = 0. This implies that in the alpha 
region one one has that for e = 0, Xhet = 0{N). If we also assume that fth{x, y) behaves quadratically for large 
negative x, we get xt/i(</>, e, A^) = N^'^ fthicliN^'^ .^N^'"^) - N^^^cji^gthieN^^^); which implies a N'^/'^ scaling of 
Xth for e = and a behavior as e~^^^ for small positive e and N — > oo. A linear behavior fth{x,y) = xhtfi{y) 
would give Xth ^ N^^^ for e = and Xth ^ for small e and N — > oo. 

Away from criticality, a simple dynamical scenario can be conjectured for the alpha relaxation in the barrier 
dominated region e < and N ^ oo. We make a crucial assumption that in a first approximation in the region 
e < 0, the reparameterization invariant fluctuations in the alpha region are equivalent to the ones of a simple 
barrier jumping process, as described for example by the Langevin equation 

^ = -iy'(g|50) + 77 (79) 
as 

in the weak noise limit 77 ^ 1. The waiting time before a jump is a random variable whose typical values are 
much larger then the time required to pass from Cp to when the jump occurs. This leads to a bistable behavior 
that naturally gives rise to 0(1) fluctuations of the overlap. At any given instant of time one has just a small 
probability of finding C{t) to be different from Cp or 0. Neglecting this probability we find a general form of 
the total susceptibility Xii9)j which is independent of the barrier. Suppose now to fix time in a way that for 
each value of the system size Cav = [(C)l- The histogram of C for different trajectories and initial conditions is 
approximately given by 

P{C\Cav) = ^S{C-Cp) + il~^)6iC) (80) 
The direct computation of XiiCav) using ((80)) leads to the simple expression 

X4Cav)/N ^Cav{Cp-Cav). (81) 

This form is independent of the barrier and valid in the whole region T < Td- As we will see in the next section, 
these hypotheses provide a good (but still not perfect) description of numerical results even at e = 0. Notice that 
in deriving (l5Tt we have just used bistability and not the detailed Langevin equation ([79]) . 

One can then try to use similar ideas to describe the alpha relaxation in finite dimensional systems. For 
example, numerical simulations in [40] suggest that bistability as described by (|80l) approximately holds for the 
probability of the local overlap qx{t) coarse grained on scales smaller than the correlation length. It would be 
then tempting to write an equation analogous to (|79|) , with an additional Laplacian term to describe the decay of 
the overlap. Unfortunately the phenomenology of such an equation |35| . which describes the competition between 
different phases in presence of disorder, is rather distant from the one of supercooled liquids. In particular such 
an equation would predict nucleation events as they could be expected in liquids, but also the fast growth of 
supercritical nuclei of low overlap that cannot be expected. 



7 Simulations 

In order to test the theoretical predictions derived above, we have simulated the 3-spin model on a random 
regular graph, i.e. a random graph of fixed connectivity z. This model involves N Ising spins Si = ±1 interacting 
through the Hamiltonian 

M 

= - ^ Jf^Sii^Sii^S^^ (82) 

where the indexes i^, ^2 j *3 ^-^e chosen at random in such a way that each spin participates exactly to z = 3AI/N 
interactions and the couplings are independent random variables taking values ±1, with P( = 1) = r. The 
properties of the model in thermodynamic limit are well known 36 . In particular if z > 4 the model exhibit a MC 
like dynamical phase transition at temperature Td and a Kauzmann ideal glass transition at a lower temperature 
Tk', both transition temperatures depend on z, but not on r. In fact it is well known that the thermodynamic 



16 



and dynamic properties of the model are independent of the choice of r for temperatures T > T^. In this range 
of temperature the model is annealed: that is, denoting Zj the partition function for a given disorder realization 
J, one has that limjv->.co -^lElog^j — limN^ca "^logEZj. This property is true in particular for the symmetric 
model with r = 1/2 and for the gauged model defined by the Nishimori condition r = (1 — e~^^)~^ l3.T^- Although 
the two versions of the model (the symmetric and the gauged one) have the same thermodynamical behavior for 
T > Tk, they may show very different finite size effects. In the following we study both versions. 

We choose to simulate the case z = 8. The analytic solution to the model |36j predicts, in the thermodynamic 
limit, a dynamical temperature Td — 1.3420(5) and the plateau value at Td equal to Cp — 0.750(5). All the 
results shown in the following have been obtained at the dynamical critical temperature T^. 

7.1 Simulations in the beta regime 

We start by showing results obtained with the symmetric model (r = 1/2). We have simulated at the dynamical 
critical temperature Td systems of size ranging from N — 60 to N — 150, with a number of samples Ns such 
that N X Ns = 1.8 • 10^. For each sample, we have obtained 2 independent equilibrium configurations (with the 
use of the parallel tempering algorithm), and, from each of these, we have evolved 2 independent trajectories 
with different thermal noises: this allow us to compute all the three different fiuctuations Xdis, Xth and Xhet (see 
below the detailed explanation for the gauged model). 




Figure 2: Three different susceptibilities Xdis, Xth and Xhet (from bottom to top) versus average correlation Cav 
for the 3-spin symmetric model. The xth data have been divided by a factor 3 in order to avoid data overlap. 
The vertical line marks the analytical value of Cp = 0.75. System sizes are N = 60, 90, 120, 150 (from bottom to 
top for each susceptibility). 

In Fig. [3] we plot these three susceptibilities Xdis, Xth and Xhet measured in the symmetric model (r = 1/2) 
for sizes N = 60,90,120,150. In the symmetric model, we can not study larger sizes due to the very large 
thermalization times. Nonetheless, even for these relatively small sizes, we can clearly see the very different size 
dependence of the susceptibilities at the critical point, Cav = Cp. As expected from the discussion in section l4Tl 
Xdis is practically size independent within error bars, while xth and Xhet are well compatible with the predicted 
scaling laws, iV^/'* and N^^"^ respectively. Unfortunately data are plagued by severe finite size effects even in the 
beta region, Cav > Cp. 

In order to reach a clearer conclusion about the scaling laws in the beta regime and at the critical point, we need 
to study larger systems and this is the reason for using the gauged model, where an equilibrium configuration 
can be generated without the long thermalization process Indeed on the Nishimori line, for a given 

interaction graph, one can first choose an arbitrary spin configuration 5'°, and then fix the coupling such that 
5° is an equilibrium configuration. This can be done by choosing the couplings as independent random variables 
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taking values ±1 with the following probability: 

exp [—l3Ja iSRi iSBj S^ti ) 
"•'■'"I*"' 2cosl.(w' ' ■ 

Since can be arbitrary, it is customary to set 5^* = 1 for all i, convention that we will adopt in the following. 
As a drawback of the method we note that since we generate at the same time initial configuration and quenched 
couplings, we cannot disentangle the contributions of these two sources of noise in the fluctuations. With this 
method we will therefore be able to compute Xth and the sum of Xhet + Xdis but not each term separately. Since 
we expect the contribution due to the quenched disorder (xdis) to be small, we will improperly call Xhet the 
latter sum. Notice that besides the choice of the random coupling J^, an additional source of quenched disorder 
comes from the choice of the random graph. It is known however that local properties of random regular graphs 
are self-averaging and consequently the effect of topology fluctuations are even smaller. 

We have simulated at the dynamical critical temperature Td systems of size ranging from N = 150 to 
N — 2400, with a number of samples Ns such that N x Ns = 3.7 • 10^. In order to be able to measure both xth 
and Xhet, for each sample a — 1, Ns we have simulated two independent trajectories s — 1,2 starting from the 
same equilibrium configuration and evolving with different thermal noises. We then for each sample and each 
trajectory measure the correlation functions Ca^sit) = 5'"(0)S'"'*(t), where the initial state is equal for the 

two trajectories. The suceptibilities are then estimated as 

Xth = ^E(^"")'-]^E^""^"" (84) 

a,s a,s 



Xhet = • (85) 




Figure 3: Thermal susceptibility Xth (lower data) and heterogeneity susceptibility Xhet (upper data) versus 
average correlation Cav for the 3-spin gauged model. The vertical line marks the analytical value of Cp. 

These susceptibilities are shown in fig. |3l The very large number of samples simulated allow to reduce the 
statistical noise and to work on very clean data. We notice that for large values of Cav the data converge, as 
expected, to a finite value in the N ^ oo limit: in this regime, dynamical fiuctuations due to heterogeneities are 
roughly one order of magnitude larger than those due to thermal noise. 

We now check finite size scaling of Xth and Xhet, which for e = read 

Xth{cl),0,N) ^ N^/^ htum^'^) , (86) 

Xhet{<t>,0,N) = N^/^ hhet{4>N^'^) . (87) 
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Figure 4: Rescaled heterogeneity susceptibility Xhet (left) and thermal susceptibility xth (right) in the gauged 
3-spin model discussed in the text at the dynamical critical temperature Td as a function of average correlation 
Cav for several sizes. The crossing point moves towards Cav — Cp (marked by a vertical line) for large N. 



In figures m we plot the rescaled susceptibilities, XthN~^^^ and XhetN"^/^ , as a function oi cf) = Cav — Cp for 
various values of N. We see that in agreement with the analytical predictions the different curves cross very close 
to C = Cp = 0.75. A more detailed analysis confirms that the crossing point tends to Cp for large N. 

Before testing the (j)N^^^^ scaling, we want to discuss the main source of finite size effects in the gauged model. 
By choosing the coupling signs independently, we have that energy fiuctuations in the initial configuration are 
0{N^^/'^), but with a rather large coefficient if compared to the thermalized symmetric model. The main 
consequence is that, even for N = 0(10'^), the vast majority of samples have either an initial energy much larger 
than Ed = —z/3ta,nh{l/Td) (and thus decorrelate very fast), either much smaller than Ed (and thus are stuck at 
C > Cp). The final effect is to have larger fiuctuations and larger finite size effects, with respect to a model where 
the initial energy is more concentrated around Ed- Given that, in the thermodynamical limit, the energy must 
converge to Ed, we introduce a new fixed- energy model where the initial energy is fixed to Ed- This is achieved 
generating samples with a fixed number of negative coupling equal to M(l — tanh l/r£j)/2. 
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Figure 5: Scaled heterogeneity susceptibility Xhet]^~^^^ in the fixed-energy 3-spin model. 
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Figure 6: Scaled thermal susceptibility Xhet -energy 3-spin model. 



We have simulated the fixed-energy model for sizes ranging from N = 300 to iV = 2400, with 2 thermal 
histories per sample and a number of samples Ns such that N x Ng — 3 • 10^. In the main panels of Figs. [S] 
and [S] we show the rescaled susceptibilities XthN~^^^ and XhetN^^^^ as a function of the average correlation for 
several system sizes. A comparison to the gauged model (see Fig. 2]) shows that finite size effects are reduced in 
this new model: indeed in figures [5] and |6] data cross exactly at Cav — Cp with almost no finite size corrections. 

We test then the (pN^^^'^ scaling in the insets of Figs. [Slandini where the same rescaled susceptibilities, are 
plotted now as a function of the scaling variable N~^/^{Cav — Cp). We note that a good scaling is observed in a 
relatively wide region around the origin. 
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Figure 7: Plot of the susceptibility ratio p = Xhet/Xth ^ function of Cav The curves cross in Cp, where, as 
seen in the inset, the slope scales as TV^/'*. 
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In figure [7] we show, as a function of Cav: the p ratio defined in Eg. ([5^ . which, for A; = 0, is 

Xhet 

As expected the curves cross at Cp. The Umit hmc-i-Cp liinAr_>.oo p{4'i^) taken in the specified order, is a direct 
measure of the variance of the random temperature entering the field theory, — (to2 + TO3). Unfortunately the 
values of N we can simulate do not allow an estimate of this limit. However we can see the scaling of slope dp/dC 
in Cp which behaves as A^^/"*, which confirms the formation of a discontinuity. 

7.2 Simulations in the alpha regime 

In order to study fluctuations in the alpha regime, the use of the flxed-energy model is mandatory. Indeed the 
finite size effects in the gauged model described above (strong fluctuations in the initial conflguration energy) 
amplify in the alpha regime, making the gauged model almost useless. The fixed-energy model, on the contrary, 
produces data showing a much better scaling. 





Figure 8: Rescaled N ^Xhet (left) and TV '^/'^xth (right) in the alpha regime. 



We note first of all that, the different scaling between fluctuations due to the heterogeneities and those due 
to the thermal noise persists in the alpha regime. Indeed Xhet and xth show different size dependence even for 
Cav < Cp, as can be seen in Fig. |S1 While the scaling of Xhet is clearly 0{N), that of xth has a smaller power: 
the data shown in Fig. [8] (right) suggest the value 3/4 argued for in section [6?2] to be an upper bound to the right 
power. 

In Fig.[S]we report data for the rescaled total susceptibility Xi/^ as a function of the average correlation. The 
qualitative behaviour seems to reflect quite well the bistable behaviour discussed in Section [51 and represented in 
Fig.|n]with the parabola Cav{Cp — Cav)- However some deviations from this behavior are expected, especially for 
Cav close to Cp, because the time to relax from Cp to is comparable to the time for entering the alpha regime 
in the fastest samples: as explained, the flnal result should be that Xi/-^ is quadratic around Cp. 

Bistability can be checked directly by looking at the histogram of correlations. In flg. [10] we report such 
histograms measured at times such that the mean correlation satisfles Cav = Cp/2; this time grows with system 
size. The shape of the histograms in Fig.[TU]is clearly made of two well separated peaks. However one may notice 
that there is a small but non-zero probability (apparently not vanishing when N — >■ 00) of finding intermediate 
values of correlation, and this is another plausible explanation for the deviations of Xi/-^ from the predicted 
parabola Cav{Cp - Cav)- 



8 Concluding remarks 

In this paper we emphasize the importance of describing fiuctuations in a reparameterization invariant form in 
glassy systems. We provide a universal theory of these fluctuations in the beta regime close to the mode coupling 
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Figure 9: Plot of the total rescaled susceptibility Xi/^ ^ ^ function of Cav in the alpha regime. The curves for 
different values of N are compared with the parabola Cav{Cp — Cav), derived under the hypothesis of a perfect 
bistability relaxation process. 
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Figure 10: Histogram of the correlation for Cav = Cp/2 on a logarithmic scale. While the data clearly show a 
bistable behavior, still a non vanishing part of the distribution between the two peaks seems to persists for large 
N. 



transition. 

There are three main physical ingredients in our theory: 
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1. Time scale separation, leading to quasi-equilibrium sampling of metastablc states. 

2. The vicinity to a dynamical critical point. 

3. Neglection of all possible non-perturbative effects. 

The first property allows us to study fluctuations through the use of constrained equilibrium measures and 
their associated glassy effective potential. These measures depend on a reference configuration which is itself 
randomly chosen with canonical distribution. In this paper we have extended the theory of the effective potential 
to study fluctuation with respect to this source of noise. We have considered the effective potential for fixed initial 
configuration as a random functional whose probability distribution relates to the one of overlap fluctuations and 
therefore to dynamical heterogeneities. In our description time is eliminated and we use the average correlation 
function as a clock. 

The second property allows us to invoke universality and to use the general form of replica field theory that can 
be obtained by symmetry considerations as an expansion around a Mode Coupling Transition point. The analysis 
of this theory leads to the identification of the relevant fluctuation modes. The effective field theory that describes 
them reduces, through tremendous simplifications to a scalar cubic field theory with a local random field term. 
The random field term is the expression of heterogeneity in the initial condition and acts as a source of disorder 
that influences the subsequent dynamics. A remarkable consequence of our description is that fluctuations with 
respect to different sources of noise show different singular behavior as the dynamic transition is approached. 

The third point is quite delicate. In deriving the equivalence with the RFIM all barrier jumping processes 
are neglected. Quite naturally one could hypothesize validity of the RFIM description beyond perturbation 
theory. However, the dynamics in the non-perturbative region, even its reparametrization invariant part, could 
be very different in the RFIM and in supercooled liquids. The analogy between liquid dynamics and the decay of 
metastable phases decay is only partially valid. In ordinary first order transition kinetics, competition between 
interface and volume free-energy leads to fast growth of supercritical nuclei. This fast growth should not be 
present in supercooled liquids. 

We tested our scenario in the favourable case of finite size scaling in mean-field models. More work will be 
needed to test the scenario in liquid models. In low dimension the critical -or pseudocritical- properties of the 
spinodal point of RFIM can hardly be computed analytically. One should therefore compare numerical results 
of simulations of liquids with numerical results on the RFIM, however the situation might be complicated by the 
fact that non-perturbative effects could be different in the two kind of systems. 

An important prediction of our theory is that the different components of the dynamical fluctuations have 
different scaling properties. Within a gaussian approximation we have found that the heterogeneous susceptibility 
is proportional to the square of the thermal one. It is not clear to us if this simple quadratic relation holds beyond 
the gaussian approximation or two not-simply related exponents describe the corresponding singularities. This 
is an important point that will need to be clarifled through numerical simulations. 
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